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(Abstract) The objective of this work is to investigate the effect of using FRP strengthening techniques to rehabilitate cracked 
structural members like beams under numerous loading conditions. Liner Elastic Fracture Mechanics was adopted in this research 
program in order to estimate Mode I Stress Intensity Factor (SIF). Standard FEM (ANSYS) was used to simulate the cracked beams 
by using 8-node quadrilateral element. The current analysis showed that SIF of the cracked structural member decreases by using FRP 
laminate because of crack bridging. In addition length, thickness and fiber volume fraction of the FRP has great influence on 
decreasing the estimated SIF. 
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1. INTRODUCTION 

Many steel structural members are in need for rehabilitation 
because they may suffer fatigue cracking or corrosion damage. 
Strengthening may also be required due to increase in their 
weights, and there is a demand to increase the flexural capacity 
of these members. The repair of the structural members by FRP 
rehabilitation techniques is one of the promising retrofit and 
common techniques. This is metallic due to its ease of 
application and the weight added to a structure is minimal. The 
cracks could propagate and cause complete failure of the 
member. The use of Fiber Reinforced Polymer (FRP) sheets to 
repair the cracked members could be a quite encouraging 
technique because replacing the whole structure could be very 
costly. This may contribute labor savings and reducing 
requirements for staging and lifting materials. Traditional 
procedures for strengthening and rehabilitation of damaged 
steel structures often require welding or bolting of the structure. 
In contrast the fatigue performance of bonded composite of 
cracked steel members has been shown to be effective up to 20 
million cycles [1]. It has also been estimated that bonded 
carbon fiber reinforced polymers CFRP strips have been used 
over 10000 times for the repair of corrosion and fatigue damage 
of aircrafts [2]. These effects illustrate the increasing 
acceptance of such technique in applications where safety and 
durability are the criteria. The goal of this research is to restore 
the strength of the cracked members, and this was investigated 
numerically by FEM. A parametric study is intended to 
demonstrate the effectiveness of a proposed repair technique 
for cracked structural members by using FRP batches. Using 



standard Finite Element (FE) program ANSYS [3], has 
simulated the cracking in the structural, members and the 
element used in the analysis is 8-node quadrilateral element. 
Three types of beams were investigated under various loading 
conditions and the numerically estimated Stress Intensity 
Factor (SIF) was calculated by using Liner Elastic Fracture 
Mechanics (LEFM) technology in order to neglect the crack tip 
plasticity [4, 5]. 

2. BACKGROUND 

Linear elastic fracture mechanics technology is based on an 
analytical procedure that relates the stress field magnitude and 
distribution in the vicinity of a crack tip to the nominal stress 
applied to the structural member, size, shape and orientation of 
the crack and material properties [4]. The fundamental 
principle of the fracture mechanics is that the stress field a 
head of a sharp crack in structural member can be represented 
by a single parameter K, (SIF) that has a unit of MN.m" 5 ' 2 . 
This parameter relates both the nominal stress level a in the 
member and the size of the crack present a .Paris and Sih [4] 
showed that in order to establish methods of stress analysis for 
cracks in elastic solid, it is convenient to define three types of 
relative movement of two crack surfaces. These displacement 
modes represent the local deformation in an infinitesimal 
element containing a crack front which is the opening mode, 
mode I, the sliding or shear mode, mode II and the tearing 
mode, mode III .It is clear that each of these modes of 
deformation corresponds to a basic type of stress field in the 
vicinity of crack tips. In any case the deformations at the crack 
tip can be treated as one or combined of these local 
displacement modes. Moreover the stress field at the crack tip 
can be treated as one or a combination of the three basic types 



of stress fields. Most practical design situations and failures 
correspond to mode I, which is the mode, considered in this 
study. 

Westergaared and Irwin [4] developed a method to show that 
the stress and displacement field in the vicinity of crack tip for a 
linear, elastic, homogeneous, and isotropic solid subjected to 
the three modes of deformation are given by: 
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where, 

IQ is stress intensity factor for mode I . 
G = Shear modulus = E/2(l+u) 

k=(3-u)/(l+p) for plane stress 
k=3-4u for plane strain 
u= Poisson's ratio 
(7 Z =0 for plane stress 
(7 z =u((7 x -(7y) for plane strain 



It may be stated that most brittle fractures occur under 
conditions of mode I loading [6], therefore the stress of 
primary intent and in most practical applications is cj y .The 
stress becomes infinite and the SIF is then a measure of the 
stress singularity at the crack tip. For u y to be a maximum in 
Eq.2, let 9=0°, therefore: 
K, 



cr =- 



■ => K I =a Y 42m 
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It is clear that at location far from crack tip, cr y decreases, and 
Ki remains constant and describes the intensity of the stress 
field just ahead of sharp crack. 

3. CASES STUDIED 

A simply supported beam of length=400mm, width=40mm 
and thickness =40mm, as shown in Figure 1. The beam is 
subjected to shear force of 5.5kN which represent 20% of the 
yield stress of the structural steel AISI 1020, is the case study 
one. While the same geometry of the beam under a tensile 
load, as shown in Figure 2, is the case study two. It is 
assumed that the beams were made of steel with Young's 
modulus of E=200* 10 3 MPa and Poisson's ratio of u=0.3. The 
investigation includes the study of SIF for both 
unstrengthened and strengthened cracked beam for the two 
cases. The crack depth of a=4mm, 8mm, 12mm, 16mm and 
20mm were studied. These values are corresponding to crack 



length to width ratio a/W of 0.1, 0.2, 0.3, 0.4 and 0.5 
respectively. The effect of stiffener dimensions on SIF, 
thickness t (20mm, 40mm and 80mm) and length L (8mm, 
20mm) were considered Fig. 3 and 4. In addition a comparison 
is made to investigate the effect of using different types of the 
fiber reinforcement (i.e. fiber glass, kevlar and carbon fiber) of 
different modulus of elasticity. This is to explore the effect of 
stiffener stiffness. 



4. FINITE ELEMENT MODELING 

The numerical modeling for all cases was based on the finite 
element method FEM by using two dimensional plain strain 
8-node quadrilateral element (Solid 86) [3] in ANSYS and is 
used to obtain well established singular stress field by shifting 
the mid-nodes one-quarter away from crack tip to compute SIF. 
Because of symmetry, only one half of the beam was modeled 
and all models were meshed with the same degree of 
refinement at the crack edge to obtain consistent result. The 
finite element boundary condition of the cracked beam is 
shown in Figure 5 and 6, while the boundary conditions for the 
strengthened cracked beam are shown in Figure 7 and 8. The 
typical FE mesh adopted in the analysis is shown in Figure 9. 



5. STRENGTHENING MATERIAL 

Three different types of fiber reinforced epoxy composite 
stiffeners are used in the study. These are fiber glass, kevlar 
and carbon fibers. Table 1 [7] show a comparison for the 
modulus of elasticity of three types. The modulus of elasticity 
and the Poisson's ratio of the used bonding adhesive (Epoxy) 
is 3.5GPa p=0.25 respectively. This is according to the data 
sheet provided by the manufacturers [8]. Because of the 
complexity of the problem investigated, the FRP batch was 
assumed to be isotropic for each of the materials. 
The composite modulus of elasticity computed by using the 
following relation [6]: 

E c =E f V f fE m V m (7) 
where E c , E f and E m are the modulus of elasticity of composite 
batch, fiber and matrix respectively, and V f and V m are the 
volume fraction of the fiber and matrix respectively. The 
volume fraction of the fibers used 25%, 50% and 75%. Table 
2 shows the effect of V f on the composite stiffener modulus of 
elasticity calculated using Eq.7. 

6. ESTIMATING OF SIF BY FEM 

Fracture mechanics generally is based on two types of 
analysis [10], namely residual strength analysis to determine 
the maximum crack size that can be tolerated, and fatigue 
crack growth analysis to calculate the time for crack growth 
from a certain initial crack size until the maximum tolerable 
crack size in order to determine the safe life. These analyses 
are usually based on the SIF intensity factor. Therefore, 
knowing the SIF is important for the fracture mechanics. 
Stress intensity factor can always be expressed as: 



K = <jy[mY (8) 

where F is a dimensionless factor which depends upon 
geometry of the body and crack size. The problem resides in the 
determination of Y. Solutions of Y for many configurations 
have been obtained. Compilations of these solutions known as 
stress intensity factor handbook [1 1] are available. Components 
that contain cracks may be subjected to one or more different 
types of load, such as uniform tensile loads, concentrated 
tensile loads, or bending loads. In such a case the total SIF can 
be found by algebraically adding the SIF that corresponds to 
each load if it is available. This method is known as the 
superposition method [10]. It can be noted that superposition 
is the addition of SIF and not of Y. The value of Y of each case 
has to be obtained separately and SIF calculated separately and 
then added [10]. Many structural configurations are so 
complicated, that a solution may not be available in handbooks. 
Consequently the finite element method was and still a 
powerful tool to determine crack tip stress fields. By finite 
element solution a complicated engineering geometry could be 
analyzed, also it has been extensively used in engineering crack 
problems, such as the determination of SIF. 

Basically there are two different approaches to determine SIF 
by the finite element method [10]. One is the direct method in 
which SIF is calculated from stress field or from 
displacement field. The second is an indirect method in 
which SIF is calculated from other parameter, the line integral 
(7-Integral). In this study the direct method is considered. The 
displacement method involves a correlation of the finite 
element nodal point displacements with the well-known crack 
tip displacement equation. For mode I, the crack tip 
displacement Eq.4 and 5 can be written as: 



U,=^J— fi(0) 
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where Ui=u and U 2 =v 
./i=cos(#/2)[fc-l + 2sin 2 (#/2)] and 

f 2 = sin(<9/2)[/fc + l-2cos 2 (6'/2)] 

By substituting a nodal point displacement U; at some point 
(r,6) near the crack tip into Eq.9 a quantity Kj can be 
calculated as: 

fix 2GU, 
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From plots of Ki as a function of r for fixed values of 8 and a 
particular displacement component, estimates of K ( can be 
made. If the exact theoretical displacement is substituted in 
Eq.10 then value of K obtained as r approaches zero would 
be the exact value of Kj. Since the finite element 
displacements are rather inaccurate at an infinitesimal distance 
from the crack tip, this process is not useable. Instead a 
tangent extrapolating of the Kj curve is used to estimate Ki 
[12]. With suitable refinement of element size the Ki curves 
obtained from element displacements rapidly approach a 
constant slope with increasing distance 
(r) from the crack tip. 



The estimation of the SIF by the tangent extrapolating 
method could be made by plotting a tangent line to the constant 
slope portion of the curve and extending the line until an 
intersection with Ki is obtained at point, which is the estimated 
value of Ki. The most accurate estimates are obtained from Ki 
curve corresponding to the normal displacement on the crack 
surface at 6=180° [10], hence Eq.10 is therefore can be written: 

K,= 2GV M (11) 
/ 2 (180°)V r 

The stress method for determining the crack tip SIF involves 
a correlation of the finite element nodal stresses with the well- 
known crack tip stress equations. 

For mode I, the crack tip stress, Eq.l, 2 and 3 can be written 
as: 
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Nodal point stresses ay 
substituted into Eq.12 
from: 



in the vicinity of the crack tip can be 
and values of Ki can be calculated 



Kj = 



(13) 



The estimation of SIF could be made by plotting K as a 
function of r for fixed value of 9 and a particular stress 
component. If the exact theoretical stresses are substituted 
into Eq.13 then the intercept of the curve with the K! axis at 
r=0 would be the exact value of Ki. Due to the inability of the 
finite element method to represent the stress singularity 
conditions at the crack tip, the K t curve for r greater than zero 
must again be extrapolated back to r=0 The extrapolated value 
of Ki at r = is the estimated K. 

Accurate evaluation of the SIF by stress method could be 
obtained from plotting K! curve corresponding to the normal 
stress on the crack tip (a y ) at angle 9=0° [13]. Hence Eq.13 
can be written: 
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7. ESTIMATING OF SIF BY FEM 

Figure 10 shows the normalized reduction in SIF using 
different batch geometries under bending condition. The V f of 
the reinforced batch was chosen to be 25%, which is the 
common V f value for hand lay-up process. The commonly used 
Vf range is approximately 30-40%. Whereas with the higher 
quality, more sophisticated and precise processes are used in 
the aerospace industry, V f approaching 70% can be successfully 
obtained [14]. Accordingly, a range of V f from 25% up to 75% 
is adopted for the batch reinforcement specification. Table 3 
demonstrates the percentage reduction for each batch 
dimension for a crack depth a/W=(0.1, 0.2, 0.3, 0.4 and 0.5).It 



is observed that the maximum drop in SIF is 77% , whereas the 
optimum batch chosen for subsequent analysis was t=0.2W, 
L=0.5W for economic consideration. For stiffened cracked 
beam subjected to tensile load, the results are shown in Figure 
11, where normalized SIF is shown. The percentage reduction 
of SIF is shown in Table 4, where the maximum drop in SIF is 
60%. Again due to cost-effective concern, an optimum 
geometry was chosen analogous to the chose dimensions for 
bending (t=0.2W, L=0.5W). Therefore these dimensions are 
adopted for the subsequent analysis. The effect of the variation 
of Vf on SIF is represented by Figure 12, and its percentage of 
reduction of SIF is obvious in Table 5. It is concluded that 
V f =25% is reasonable for rehabilitation, although, pre-preg 
strips can be used to add more stiffness. Using different types of 
fibers, the impact on the SIF is shown in Figure 13 and the 
percentage drop in SIF is presented in Table 6. This is for the 
flexural loading case. Although carbon fiber showed the 
highest reduction, but this is due to its higher stiffness. The 
effect of using different V f and different types of fibers for 
tensile case are shown in Figure 14 and 15, which are 
explained by the Table 7 and 8. It is clear that the maximum 
reduction in SIF by was 56% for V f =75%, whereas the 
reduction was 67% for the carbon fiber stiffener of V f =25%. 

8. CONCLUSION 

The FE technique is used to estimate mode I SIF of cracked 
beam for both tensile and bending load. The following 
conclusions illustrate a summary of the results: 

1. Fiber volume fraction of 25% is the major reduction in SIF, 
although 50%-75% showed more reduction in SIF, but is not 
significant. 

2. For economical point of view V f =25% is recommended 
since the reduction in the SIF using higher Vf has small effect 
in the SIF. 

3. Glass fiber reinforced stiffener demonstrated a proper 
stiffness to be utilized for the rehabilitation of cracked beams 
in term of cost effective technique and to avoid the corrosive 
galvanic effect between carbon and the steel. 
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Table I.Modulus of elasticity for different fibers 



Modulus of Elasticity GPa 


Epoxy 


E-Glass 


Kevlar 


Steel 


Carbon 


3.5 


76 


124 


205 


640 



Table 2. FRP's Modulus of Elasticity versus Vf% 



Modulus of Elasticity GPa 


V f % 


E-Glass 


Kevlar 


Carbon 


25 


21.6 


33.6 


162.6 


50 


39.75 


63.7 


321.7 


75 


57.87 


93.87 


480 



Table 3.Percentage reduction of SIF using different batch 

GEOMETRY, VF=25%, BENDING 



a/W 


t=0.2W 
L=0.5W 


t=0.2W 
L=W 


t=0.2W 
L=2W 


t=0.5W 
L=W 


t=0.5W 
L=2W 


0.1-0.5 


24-65% 


28-69% 


29-70% 


31-71% 


41-77% 



Table 4.Percentage reduction of SIF using different batch 

GEOMETRY, VF=25%, TENSION 



a/W 


t=0.2W 
L=0.5W 


t=0.2W 
L=W 


t=0.2W 
L=2W 


t=0.5W 
L=W 


t=0.5W 
L=2W 


0.1-0.5 


30-50% 


31-51% 


32-52% 


34-54% 


38%-60% 



Table 5.Percentage reduction in SIF using different Vf 
ratio(Glass fiber), Bending case 



a/W 


V t =25% 


Vf=50% 


V t =75% 


0.1-0.5 


24-65% 


30-68% 


34-70% 



table 6.percentage reduction of sif by using different types of 
fibers(Vf=25%), Bending case 



a/W 


Glass 
fiber 


Kevlar 


Carbon 


0.1-0.5 


24-65% 


29-69% 


50-77% 



Table 7.Percentage reduction in SIF using different Vf 
ratio(Glass fiber), Tensile case 



a/W 


V t =25% 


Vf=50% 


V t =75% 


0.1-0.5 


30-50% 


34-53% 


38-56% 



Table 8.Percentage reduction of SIF by using different types of 
fibers(Vf=25%), Tensile case 



a/W 


Glass fiber 


Kevlar 


Carbon 


0.1-0.5 


30-50% 


35-53% 


53-67% 
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Figure 1. Cracked beam subjected to bending load. 




Figure 2. Cracked beam subjected to tensile force. 
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Figure 3. Cracked strengthen beam subjected to bending load. 
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Figure 4. Cracked strengthened beam subjected to tensile load. 
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Figure 6. FE simulation of the half length of a cracked beam 
subjected to tensile load. 




Figure 7. FE simulation of the half length of cracked strengthened 
beam subjected to bending force. 




Figure 8. FE simulation of the half length of a cracked 
strengthened beam subjected to tensile load. 
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Figure 9. FE mesh of the half length of a cracked strengthened 
beam. 



Figure 5. FE simulation of the half length of cracked beam 
subjected to bending force. 



-6- 




Figure 14.Effect of V f % on SIR 
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Figure 15. Variation of SIF by fibers type (V f =25%). 



